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1.  Introduction 


The  Gilbert-Johnson-Keerthi^  (GJK)  algorithm  is  a  very  effieient  algorithm  for  de¬ 
termining  the  minimum  distance  between  2  convex  objects.  GJK  also  doubles  as  a 
very  fast  collision  detection,  which  has  proven  extremely  useful  in  modeling  work.^ 
In  particular,  while  modeling  the  system  effectiveness  of  different  threats,  sensors, 
and  armor  packages,  we  required  a  way  of  very  quickly  determining  if  a  threat  and 
counter-munition  had  collided,  if  a  threat  had  impacted  a  vehicle,  or  if  multiple 
counter-munitions  had  collided,  and  so  on.  All  of  these  problems  required  a  very 
fast  algorithm  to  determine  if  2  objects  intersected  where  the  objects  may  not  be 
simple  lines,  cubes,  cylinders,  or  other  polygons. 

While  GJK  was  useful  in  our  work,  implementing  GJK  proved  to  be  extremely  diffi¬ 
cult.  The  original  GJK  paper^  takes  an  algebraically  intensive  approach  to  compute 
results,  which  becomes  difficult  to  implement.  Alternatively,  GJK  can  be  imple¬ 
mented  using  more  geometrical  concepts,  which  creates  a  more  intuitive  and  under¬ 
standable  algorithm.  Numerous  tutorials  have  been  posted^"^  adopting  the  geomet¬ 
rical  approach.  However,  we  encountered  the  following  difficulties: 

•  The  tutorials  do  not  agree  with  each  other.  The  basic  concepts  are  the  same 
across  various  tutorials,  but  the  implementation  details  conflict  with  each 
other.  Certain  cases  may  be  ignored  or  handled  inconsistently. 

•  Most  tutorials  only  focus  on  2  dimensions.  The  3-dimensional  (3-D)  case, 
however,  is  not  a  straightforward  generalization  of  the  2-dimensional  (2-D) 
case. 

•  Very  little  is  written  about  the  assumptions  involved  in  the  geometrical  case. 
This  leads  to  difficulties  since  certain  portions  of  the  algorithm  require  very 
specific  orderings. 

These  problems  may  not  be  recognized  by  the  tutorial  authors  themselves  since 
GJK  is  an  iterative  algorithm  that  improves  results  from  one  iteration  to  the  next 
until  an  answer  is  found.  The  iterative  nature  of  GJK  can  mask  many  bugs  in  an 
implementation  because  the  implementation  can  self-correct  in  later  iterations.  That 
is,  if  a  bug  exists  in  a  particular  case,  the  implementation  may  cause  GJK  to  search 
for  results  in  a  wrong  direction.  However,  later  iterations  could  eventually  bring  the 
implementation  back  to  the  correct  direction.  Thus,  an  implementation  can  have  a 
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bug,  but  still  pass  many  test  cases,  resulting  in  the  bug  staying  in  the  code  for  a  long 
time. 

These  problems  have  led  us  to  record  the  details  of  our  implementation  of  the  GJK 
algorithm.  We  discuss  the  assumptions  going  into  each  portion  of  the  implementa¬ 
tion,  as  well  as  why  certain  orderings  are  required.  Given  the  difficulties  above,  we 
cannot  state  that  our  implementation  is  bug  free.  However,  we  do  discuss  our  efforts 
at  verifying  our  implementation  against  an  outside  source. 

2.  Mathematical  Preliminaries 

We  first  discuss  a  few  mathematical  concepts  that  are  needed  later.  Readers  familiar 
with  this  material  can  simply  pass  to  the  next  section. 

We  use  the  notation  that  A  refers  to  a  point  in  space,  while  A  refers  to  the  vector 
from  the  origin  to  the  point  A.  AB  =  B  —  A  =  vector  from  point  A  to  point  B.  We 
use  Ai  to  refer  to  the  f-th  component  of  vector  A.  We  say  that  A  is  a  zero  vector  if 
all  A,  =  0. 

Points  and  vectors  may  be  gathered  into  sets,  which  we  indicate  with  script  letters, 
e.g.,  si ,  Two  points  can  also  form  a  line  segment,  which  we  indicate  as  AB. 
Three  (non-colinear)  points  form  a  triangle,  which  we  indicate  as  /AABC. 

2.1  Dot  and  Cross  Products 

The  (scalar)  dot  and  (vector)  cross  product  are  the  2  most  basic  product  operations 
involving  vectors.  The  dot  product  in  Cartesian  coordinates  is 

\B  =  '£A,Bi,  (1) 

i 

Geometrically,  A  •  B  >  0  means  that  A  has  a  component  in  the  same  direction  as  B 
(and  vice-versa).  Informally,  we  can  say  that  the  2  vector  point  roughly  in  the  same 
direction  if  the  dot  product  is  positive.  A  ■  B  <  0  means  the  vectors  point  in  opposite 
directions.  A  •  B  =  0  means  the  2  vectors  are  perpendicular  provided  that  A  and  B 
are  not  zero  vectors. 

The  cross  product  is  vector  whose  z-th  component  is  defined  as 


(A  X  B);'  —  Eijj^AjBi^., 


(2) 
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where  is  the  Levi-Civita  symbol  whose  value  is  1  if  i jk  is  an  even  permutation 
(i.e.,  123,  231,  or  312),  —1  if  z  jk  is  an  odd  permutation  (i.e.,  321,  213,  or  132),  and 
0  otherwise.  The  x  component  of  the  cross  product  is  thus  AyB^  —  A  Jiy.  Geometri¬ 
cally,  A  X  B  is  a  vector  that  is  perpendicular  to  both  A  and  B.  This  interpretation  is 
very  important  in  Cartesian  coordinates  where  2  linearly  independent  vectors  define 
a  plane.  In  this  case,  the  cross  product  of  the  vectors  would  define  a  vector  normal 
to  the  plane. 

For  the  dot  product,  the  order  of  the  vectors  is  unimportant  since  A  •  B  =  B  ■  A. 
This  is  not  true  for  the  cross  product  where  changing  the  order  inverts  the  sign,  i.e., 
A  X  B  =  — B  X  A.  Geometrically,  this  implies  that  a  plane  has  2  unit  normals  that 
point  in  opposite  directions.  This  is  very  important  for  the  GJK  algorithm  where 
we  take  the  cross  products  in  very  particular  orders  so  that  the  normals  point  in  the 
direction  we  require. 

In  some  instances,  we  require  the  vector  cross  product  of  3  vectors,  A,B,  and  C. 
Parentheses  may  be  used  to  specify  which  cross  product  should  be  performed  first. 
If  no  parentheses  are  given,  then  the  cross  products  proceed  from  left  to  right.  That 
is,  A  X  B  X  C  is  evaluated  as  (A  x  B)  x  C. 

2.2  Convex  Sets  and  Convex  Hulls 

GJK  works  on  convex  objects  (and  sets).  Let  i?/  be  a  set  and  let  Pi  and  P2  be  any 
2  points  in  £/.  If  all  the  points  on  the  shortest  path  (a  line  in  Euclidean  geometry) 
between  Pi  and  P2  are  also  in  sA ,  then  .<2/  is  a  convex  set. 

Figure  1  shows  some  examples  of  convex  sets.  A  circle  is  not  convex  because  it 
does  not  include  any  interior  points.  Thus,  any  chord  of  the  circle  is  not  included  in 
the  set,  which  means  the  set  cannot  be  convex.  By  contrast,  a  solid  disk  is  convex 
because  the  disk  includes  the  interior  points.  In  general,  any  hollow  object  (circle, 
empty  cube,  etc.)  is  not  convex  because  the  interior  points  are  not  included.  Only 
solid,  “filled  in,”  objects  can  be  convex.  Note,  however,  that  not  all  solid  objects  are 
convex  as  the  crescent  in  Fig.  1  demonstrates. 

GJK  only  works  with  convex  sets  and  objects.  This  restriction  means  that  GJK 
may  not  be  able  to  determine  if  2  objects  with  holes  in  them  intersect.  In  practical 
terms,  this  restriction  is  rarely  an  inconvenience  as  most  objects  can  be  regarded  as 
convex  for  the  purposes  of  determining  an  intersection.  This  is  true  for  our  modeling 
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Fig.  1  Examples  of  convex  and  non-convex  objects.  A  circle  is  not  convex  because  any  chord 
is  not  in  tbe  set.  Similarly,  a  solid  crescent  is  not  convex  because  a  line  segment  for  points  in 
tbe  wedge  area  are  not  included.  A  solid  disk  is  convex. 


work,  where  we  ehoose  eylinders  and  other  eonvex  objeets  to  represent  threats  and 
vehieles.  This  is  also  true  for  objeets  that  do  have  holes  in  them  provided  those 
objeets  interseet  with  other  objeets  larger  than  the  holes.  For  example,  eonsider  a 
net  thrown  over  a  person.  Strietly  speaking,  neither  the  net  nor  the  person  ean  be 
regarded  as  eonvex.  The  net  has  a  lot  of  empty  spaees  between  threads  while  the 
person  could  raise  her  limbs  out  from  the  body,  creating  empty  spaces.  However, 
for  the  purposes  of  throwing  a  net  over  the  person,  we  can  regard  both  as  being 
convex  provided  the  spaces  in  the  net  are  too  small  to  pass  over  the  person.  In  such 
a  scenario,  the  net  will  still  drape  over  the  person  as  if  both  were  convex. 

We  do  not  necessarily  need  to  keep  track  of  all  the  points  in  the  convex  set  . 
Since  any  line  segment  joining  2  points  in  si  must  also  be  in  si ,  there  must  be 
some  minimal  number  of  points  on  the  outer  edge  of  si  such  that  we  can  draw  line 
segments  between  those  points  and  still  recover  all  of  si .  That  is,  there  is  some  set 
si'  C  si  such  that  si'  is  convex  and  contains  the  region  bounded  by  si .  The  points 
in  si'  form  the  outer  envelope  (or  hull)  of  si .  Since  si'  is,  by  construction,  convex, 
we  term  this  set  the  convex  hull  of  si ,  which  we  write  as  si'  =  Com{si). 

Figure  2  shows  a  useful  visualization  of  a  convex  hull  based  on  the  rubber  band 
analogy.  Suppose  we  have  a  finite  set  of  2-D  points.  We  first  place  a  pin  at  each 
point  and  stretch  a  giant  rubber  band  around  all  the  pins.  If  we  let  go  of  the  rubber 
band,  it  will  collapse  around  the  pins  on  the  exterior  of  the  set.  The  rubber  band 
would  then  form  the  convex  hull.  This  analogy  also  demonstrates  that  for  any  finite 
set  of  points,  we  can  always  determine  a  convex  hull  covering  all  the  points. 
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Fig.  2  Sketch  of  the  ruhher  hand  analogy  for  computing  convex  hulls.  We  stretch  a  ruhher 
hand  around  a  set  of  points  and  let  the  ruhher  hand  go.  The  hand  will  snap  around  the  outer¬ 
most  points,  forming  the  convex  hull. 


For  a  solid  disk,  the  convex  hull  contains  all  the  points  on  the  circumference  of  the 
disk,  i.e.,  the  bounding  circle.  This  may  be  confusing  at  first  since  a  circle,  by  itself, 
is  not  convex.  However,  the  circle  is  the  convex  hull  of  a  disk,  which  is  convex. 
The  crucial  insight  is  to  realize  that  the  convex  hull  of  a  set  bounds  the  original 
set.  If  a  circle  is  the  convex  hull,  then  the  original  set  included  all  the  points  on  the 
circumference  of  the  circle  as  well  as  the  interior  points,  i.e.,  the  solid  disk. 

A  convex  hull  generally  simplifies  algorithms  since  we  can  analyze  a  smaller  subset. 
A  triangle,  for  example,  is  a  convex  object  containing  all  the  interior  points  and  the 
perimeter  connecting  the  vertices.  Any  analysis  we  perform  would  need  to  account 
for  all  those  points.  However,  we  can  specify  the  convex  hull  of  the  triangle  using 
only  the  3  vertices  of  the  triangle  since  those  3  points  bound  the  triangle  region.  By 
using  the  convex  hull,  we  only  need  to  analyze  3  points.  Algorithms,  such  as  the 
Graham  scan,  exist  to  efficiently  compute  the  convex  hull  of  any  finite  set. 

For  a  finite  set  of  N  points,  the  convex  hull  may  involve  up  to  N  points.  If  N  is  very 
large,  keeping  all  N  points  may  be  impractical.  Depending  on  the  application,  we 
may  not  need  to  keep  all  the  points.  The  theorem  of  Caratheodory  states  that  given 
a  set  £/  consisting  of  points  in  d  dimensions  and  a  point  P  G  Com{£^),  there  exists 
a  subset  C  ^  such  that  has,  at  most,  d-\-\  points  and  P  E  For 

3  dimensions,  this  means  we  need,  at  most,  4  points  from  to  construct  a  convex 
hull  containing  a  given  P. 

The  algebraic  approach  to  GJK  involves  computing  convex  hulls  of  small  sets.  The 
geometric  approach  to  GJK  does  not  directly  compute  the  convex  hulls.  Instead,  the 
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geometric  approach  uses  convex  hulls  and  the  theorem  of  Caratheodory  indirectly, 
as  is  seen  in  the  next  section. 

2.3  Minkowski  Sum  (and  Difference) 

Suppose  we  have  2  sets,  and  consisting  of  vectors.  We  can  form  a  new  set 
by  adding  each  vector  in  to  each  vector  in  SB,  and  another  new  set  by  taking 
their  difference.  This  is  the  definition  of  the  Minkowski  Sum  and  Difference  sets. 
Formally, 

Msum(.<2ii^,.^)  =  {X  + Y  I  X  G  A  Y  G  (3) 

Mdiff(.(S^,.^)  =  {X-Y|XG.<a^AYG^}  (4) 

Formally  speaking,  Msum  and  Mdiff  only  apply  to  sets  of  vectors.  However,  in 
Cartesian  geometry,  the  concepts  can  be  applied  just  as  easily  to  points  in  space. 
Thus,  we  often  use  Msum  and  Mdiff  with  respect  to  points.  Figure  3  shows  an 
example  of  the  Minkowski  Difference  for  2  triangles. 


Fig.  3  Two  example  triangles  and  their  Minkowski  Difference 


The  Minkowski  Difference  has  the  property  that  if  2  sets  intersect,  then  Mdiff  will 
contain  the  origin.  This  is  trivial  to  see  since  2  sets  intersecting  means  they  must 
have  at  least  1  point  in  common,  so  Mdiff  will  include  the  origin.  This  property  is 
important  for  GJK’s  operation. 

Rather  than  simple  sets,  we  are  concerned  with  objects,  such  as  a  triangle,  a  sphere, 
or  any  convex  polygon.  The  Minkowski  Sum  and  Difference  are  still  defined  and 
are  convex,  provided  the  objects  are  convex.  If  the  2  objects  intersect,  the  origin  is 
contained  within  the  Minkowski  Difference.  More  specifically,  the  origin  lies  some¬ 
where  within  the  volume  bounded  by  the  convex  hull  of  the  Minkowski  Difference. 
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GJK  decides  if  2  objects,  and  intersect  by  determining  if  Mdiff(^,  M)  con¬ 
tains  the  origin  or  not.  GJK  does  not  need  to  compute  the  entire  Minkowski  Differ¬ 
ence,  which  would  be  impractical  for  large  objects  or  objects  like  spheres.  Instead, 
GJK  relies  upon  support  functions  to  determine  points  on  the  convex  hull  of  the 
Minkowski  Difference. 

2.4  Support  Functions 

The  support  function  for  a  convex  object  takes  a  direction  and  returns  the  vector 
to  a  point  in  the  object  that  is  farthest  in  the  given  direction.  Support  functions  do 
not  have  to  return  a  unique  point.  Consider  triangle  shown  in  Fig.  3.  There  is 
no  vertex  point  exactly  in  the  x  direction,  which  means  the  support  function  could 
either  interpolate  between  the  points  (0,-1)  and  (0, 1),  or  simply  return  either  of 
those  vertices  as  both  are  equally  far  in  the  x  direction. 

If  the  object  consists  of  a  small  number  of  points  on  the  convex  hull,  such  as  the 
vertices  of  a  triangle,  then  the  support  function  simply  returns  a  vector  to  the  point 
whose  dot  product  with  the  direction  is  maximal.  Mathematically,  for  object  G  and 
direction  d, 

support(  ^,d)  =  P|Fg  (^AP-d>Q-d  VQ  G  ^  (5) 

Computationally,  this  is  straightforward  as  we  simply  take  the  dot  product  of  each 
point  in  the  object  with  d  and  return  the  point  that  generates  the  largest  dot  product. 

If  an  object  contains  numerous  (or  infinite)  number  of  points,  then  the  support  func¬ 
tion  needs  to  be  defined  in  a  more  appropriate  manner.  As  an  example,  consider  a 
sphere  with  radius  R  and  center  at  point  C.  The  support  function  would  then  be 

support(,5^,  d)  =  C  -f  I?  .  (6) 

In  this  case,  the  support  function  has  a  rather  simple  form.  Other  objects  may  have 
more  complex  support  functions. 

The  support  functions  for  Msum  and  Mdiff  can  become  quite  complex  as  they  in¬ 
volve  the  sum  and  difference  of  arbitrary  objects.  Fortunately,  there  are  simpler 
ways  of  computing  the  support  functions.  For  Msum,  we  note  that  the  support  func¬ 
tion  always  returns  a  point  farthest  in  the  given  direction.  Since  Msum  involves 
adding  points  from  the  2  sets,  we  can  obtain  the  farthest  point  in  Msum  by  getting 
the  point  farthest  along  the  direction  in  and  adding  the  point  farthest  along  in 
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In  other  words, 


support(Msum(^,  ^),d)  =  support(^/,  d)  +  support(e^,  d)  (7) 

where  d  is  the  direction. 

Unfortunately,  the  support  function  for  Mdiff  is  not  as  simple.  We  cannot  go  as  far 
as  possible  in  set  and  subtract  the  farthest  point  in  The  resulting  point  may 
not  be  the  farthest  point  in  Mdiff (.2/,  ^).  Instead,  we  use 

support(Mdiff(i2/,  e^),d)  =  support(.c/,  d)  —  support(e^,  — d)  (8) 

That  is,  we  go  as  far  as  we  can  along  d  in  set  and  subtract  off  the  point  in  ^  that 
is  farthest  in  the  opposite  direction,  — d. 

As  an  example,  we  consider  the  sets  shown  in  Fig.  3,  and  compute  support(Mdiff('^,  ^),x). 
support(^,  x)  returns  a  vector  to  1  of  the  points  on  the  .r  =  —  1  line.  Support(^,  — x) 
returns  a  vector  to  1  of  the  points  on  the  x  =  1  line.  For  both  cases,  multiple  points 
could  be  selected,  so  the  support  functions  for  the  sets  are  free  to  choose  any  of 
them.  Regardless,  subtracting  the  two  vectors  gives  a  vector  to  a  point  with  x  =  —2, 
which,  as  shown  in  Fig.  3,  lies  farthest  in  the  x  direction  of  Mdiff('^,  ^). 

GJK  uses  support  functions  to  choose  points  in  Mdiff  without  having  to  compute 
Mdiff  itself.  This  allows  GJK  to  be  used  for  any  convex  object  provided  the  object 
has  a  support  method. 

2.5  Simplices 

A  simplex  is  the  convex  generalization  of  a  triangle  to  an  arbitrary  number  of  di¬ 
mensions.  A  0-simplex  consists  of  a  single  point.  A  1-simplex  is  the  line  segment 
between  the  2  endpoints.  A  2-simplex  is  a  triangle  among  3  points.  A  3-simplex  is  a 
tetrahedron.  In  general,  a  J-simplex  has  d  dimensions  and  consists  of  J  -|-  1  points. 
Importantly,  a  simplex  is  always  the  convex  hull  of  the  given  number  of  vertices. 

GJK  uses  simplices  to  help  determine  if  the  origin  is  located  within  Mdiff.  GJK 
constructs  a  simplex  by  using  the  support  function  to  select  points  along  the  convex 
hull  of  Mdiff.  Since  the  resulting  simplex  has  a  simple  form  (a  point,  line  segment, 
triangle,  or  tetrahedron)  GJK  can  quickly  determine  if  the  origin  is  contained  in  the 
simplex  or  not.  Since  the  theorem  of  Caratheodory  guarantees  that  any  point  can  be 
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contained  within  a  convex  hull  of,  at  most,  4  points  (for  3  dimensions),  GJK  only 
needs  a  0-,  1-,  2-,  or  3-simplex. 

3.  GJK  Algorithm 

GJK  determines  if  2  objects,  and  e^,  intersect  (or  collide).  We  first  provide  a 
quick  overview  of  how  GJK  operates  and  then  proceed  to  a  more  in  depth  discus¬ 
sion. 

3.1  Basic  Operation  of  GJK 

3.1.1  Overview  of  GJK 

If  and  intersect,  then  their  Minkowski  Difference,  Mdiff(i2/,  will  contain 
the  origin.  Thus,  GJK  searches  Mdiff  to  determine  if  the  origin  is  contained  within 
it.  However,  Mdiff  could  be  very  complex  making  its  direct  computation  intractable. 
Instead,  GJK  makes  use  of  the  theorem  of  Caratheodory,  which  says  that  if  a  point 
P  lies  within  any  convex  hull,  then  P  also  lies  within  a  reduced  convex  hull  con¬ 
sisting  of,  at  most,  4  points  (in  3  dimensions)  of  the  larger  convex  hull.  Using  the 
theorem,  GJK  focuses  on  constructing  simplices  from  the  points  on  the  convex  hull 
of  Mdiff.  The  simplices  have  very  simple  forms,  being  either  a  point  (0-simplex),  a 
line  segment  (1-simplex),  a  triangle  (2-simplex)  or  a  tetrahedron  (3-simplex),  which 
allow  GJK  to  quickly  make  determinations  about  whether  the  origin  lies  within  the 
simplices. 

The  problem  would  be  extremely  simple  if  the  proper  simplex  were  known  at  the 
start.  Since  the  proper  simplex  is  not  known,  GJK  must  iteratively  update  the  sim¬ 
plex  to  use  additional  points  while  disregarding  points  that  are  no  longer  useful. 
GJK  determines  what  points  to  add  and  what  points  to  remove  based  on  the  direc¬ 
tion  from  the  current  simplex  to  the  origin. 

3.1.2  Two  Examples  of  GJK  Operating 

As  an  example,  we  consider  two  Mdiff  sets  shown  in  Fig.  4.  We  first  consider  the 
set  on  the  left,  which  was  made  by  two  intersecting  objects,  and  thus  includes  the 
origin.  GJK  begins  by  selecting  a  point  on  the  convex  hull  of  the  Mdiff  by  using 
the  support  function.  That  leaves  a  question  as  to  what  direction  should  be  passed 
to  the  function.  Since  GJK  is  just  beginning,  any  direction  will  be  satisfactory,  so 
we  simply  choose  the  x-direction.  The  support  function  returns  point  A,  which  GJK 
adds  to  the  simplex.  Since  A  is  the  only  point  in  the  simplex,  GJK  now  has  a  0- 
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simplex.  To  locate  the  next  point,  GJK  chooses  a  direction  that  points  from  the 
0-simplex  toward  the  origin.  The  support  function  returns  the  point  E,  which  is 
added  to  the  simplex,  creating  a  1-simplex  consisting  of  AE,  the  line  segment  from 
Ato  E.  Note  that  the  point  E  was  not  returned  by  the  support  functions  because  E 
lies  slightly  farther. 


Fig.  4  Two  example  Mdiff  sets  demonstrating  how  GJK  operates.  The  set  on  the  left  includes 
the  origin  while  the  set  on  the  right  does  not  (the  origins  are  indicated  hy  the  small  plus  signs). 
The  dashed  lines  indicate  the  segments  added  to  the  simplex  hy  GJK. 


GJK  next  chooses  a  direction  from  the  1-simplex  toward  the  origin.  The  support 
function  chooses  point  G,  which  GJK  adds  to  the  simplex,  forming  a  2-simplex  (a 
triangle).  Note  that  point  E  was  not  selected.  GJK  chose  a  direction  from  the  seg¬ 
ment  AE,  not  from  either  of  the  end  points.  The  direction  from  the  segment  toward 
the  origin  is  not  strictly  downward.  Instead,  the  direction  includes  a  horizontal  com¬ 
ponent,  which  causes  the  support  function  to  choose  G.  Section  3.2.2  shows  how  to 
compute  this  direction. 

GJK  now  has  a  2-simplex  consisting  of  the  triangle  formed  by  the  points  A,  E,  and 
G.  Since  the  origin  is  not  within  the  triangle,  GJK  has  to  select  another  point.  In  2 
dimensions,  the  2-simplex  is  the  largest  simplex  required,  so  GJK  removes  a  point 
before  adding  a  new  point.  Given  the  example,  the  origin  lies  on  the  other  side  of 
the  segment  E G,  so  point  A  is  no  longer  required.  GJK  removes  A  from  the  simplex, 
creating  a  1-simplex  consisting  of  the  points  E  and  G.  GJK  then  chooses  the  next 
point  by  moving  from  the  1-simplex  toward  the  origin  and  obtains  point  E  from  the 
support  function.  GJK  adds  E  to  the  simplex,  creating  a  new  2-simplex  consisting 
of  £■,  G,  and  E .  The  origin  is  contained  within  the  simplex,  causing  GJK  to  stop  and 
report  an  intersection. 
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Now  consider  the  Mdiff  set  shown  on  the  right-hand  side  of  Fig.  4.  This  set  does 
not  eontain  the  origin,  indieating  that  the  objeets  do  not  interseet.  GJK  proeeeds  in 
a  similar  fashion  as  above.  The  first  point  added  is  again  A.  Then  GJK  moves  from 
A  toward  the  origin  and  obtains  point  £.  GJK  now  has  a  1 -simplex  eonsisting  of  the 
ling  segment  from  A  to  E.  So  far  nothing  is  different  than  the  earlier  example.  Now, 
however,  GJK  moves  from  the  line  segment  AE  toward  the  origin,  whieh  returns 
point  E .  GJK  now  has  a  2-simplex  with  points  A,  E,  and  E .  Sinee  this  2-simplex 
does  not  eontain  the  origin,  GJK  removes  point  A  and  tries  again  to  find  another 
point.  No  further  points  ean  be  added  sinee  no  points  between  the  origin  and  the 
segment  EE  exist  in  Mdiff.  GJK  eannot  find  the  origin  within  Mdiff  and  so  reports 
no  interseetion. 

These  two  examples  demonstrate  the  operation  of  GJK  as  well  as  how  GJK  termi¬ 
nates.  GJK  operates  by  ealling  the  support  funetions  to  seleet  points  on  the  eonvex 
hull  of  Mdiff(.</,ef^).  GJK  determines  the  direetion  to  pass  to  the  support  fune- 
tion  by  eonsidering  the  eurrent  simplex  and  the  origin.  GJK  determines  the  feature 
(point,  line  segment,  or  triangle  faee)  of  the  simplex  that  lies  elosest  to  the  origin. 
Any  point  not  involved  in  that  feature  is  removed  from  the  simplex.  GJK  then  moves 
from  the  remaining  points  toward  the  origin  to  loeate  the  next  point  to  add. 

3.1.3  Termination  Conditions 

GJK  terminates  either  when  the  origin  is  loeated,  or  when  GJK  runs  out  of  points 
to  add  to  the  simplex.  The  latter  eondition  ean  be  trieky  to  implement  beeause  the 
support  funetions  are  not  required  to  return  unique  points.  If  the  set  eontains  nu¬ 
merous  points  all  equally  distant  along  some  direetion,  the  support  funetion  is  free 
to  return  any  of  them.  Thus,  we  may  not  be  able  to  determine  whether  we  have  run 
out  of  points  sinee  the  support  funetion  eould  eontinue  to  supply  new  and  different 
points. 

We  require  a  better  termination  eondition  in  the  event  that  the  origin  is  not  eontained 
within  the  simplex.  We  note  that  the  support  funetions  return  a  point  farthest  in  the 
given  direetion.  If  we  ehoose  a  direetion  that  points  toward  the  origin  and  move  as 
far  as  we  ean  in  that  direetion,  then  the  support  funetion  will  return  a  point  that  will 
be  on  the  other  side  of  the  origin  if  the  origin  is  eontained  within  Mdiff.  That  is,  we 
always  eross  over  the  origin,  if  the  origin  lies  within  Mdiff.  We  see  this  in  left  side 
of  Fig.  4,  where  we  erossed  the  origin  eaeh  time  we  added  a  new  point.  The  origin 
was  erossed  when  moving  from  A  to  E,  from  AE  to  G,  and  again  when  moving  from 
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EG  io  F .  By  contrast,  for  the  right  side  of  Fig.  4,  we  did  not  eross  the  origin  when 
moving  from  Aio  E  sinee  the  origin  was  not  within  Mdiff. 

To  see  how  we  ean  eheek  for  this  eondition,  we  examine  moving  from  point  A  to 
point  E  for  the  examples  in  Fig.  4.  A  eloseup  of  the  situation  is  shown  in  Fig.  5.  GJK 
moved  from  point  A  along  direetion  d.  The  support  funetion  in  both  eases  returned 
point  £. 


Fig.  5  Closeup  view  of  moving  from  point  A  to  point  E  for  the  examples  from  Fig.  4.  For 
both  cases,  GJK  computed  the  direction  d  to  search  for  a  new  point,  and,  through  the  support 
function,  found  point  E.  On  the  left,  the  origin  was  crossed  while  moving  from  A  to  E,  while  on 
the  right  the  origin  was  not  crossed. 


The  left  side  of  Fig.  5  shows  the  that  the  origin  was  erossed  while  moving  from 
point  A  to  point  E,  while  the  right  side  shows  the  ease  where  the  origin  was  not 
erossed.  In  both  oases,  we  ean  draw  the  veetor  E,  whioh  points  from  the  origin  to 
point  E.  From  the  figure,  we  see  that  when  the  origin  is  erossed,  d  has  a  positive 
projeotion  along  E.  On  the  other  hand,  if  the  origin  is  not  erossed,  then  d  has  a 
negative  projeetion  along  E. 

This  insight  provides  us  with  the  termination  eondition.  At  eaoh  iteration,  GJK  oom- 
putes  a  direetion,  d,  to  searoh  for  a  new  point  to  add  to  the  simplex,  d  is  passed  to 
the  support  funetion,  whieh  returns  a  new  point  P.  If 

P-d<0  (9) 

then  the  origin  was  not  erossed  when  moving  to  P  and  therefore,  the  origin  is  not 
loeated  within  Mdiff.  GJK  ean  report  no  interseetion  for  this  ease. 

Note  that  the  eondition  given  in  Eq.  9  uses  the  less-than  operator  and  not  less-than- 
or-equal.  It  is  possible  that  d  is  perpendieular  to  OP,  in  whioh  ease,  OP  ■  d  =  0.  For 
this  ease,  the  origin  oould  be  loeated  on  the  edge  of  Mdiff,  inside  of  Mdiff,  or  on 
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the  outside  of  Mdiff.  We  cannot  say  for  certain  which  case  applies,  so  GJK  cannot 
terminate  yet. 

Finally,  we  note  that  this  termination  condition  requires  a  simplex  with  at  least  1 
point  in  it  already.  In  order  for  the  concept  of  crossing  the  origin  to  make  any  sense, 
we  must  start  from  a  point  on  the  outside  of  Mdiff  and  move  to  another  point  in  the 
direction  of  the  origin.  Thus,  the  simplex  must  have  2  points  for  the  condition  to  be 
applied.  Generally,  this  is  only  an  issue  at  the  very  beginning  of  the  algorithm  when 
we  have  not  added  any  points  to  the  simplex  yet. 

3.1,4  GJK  Algorithm 

With  the  above  discussions  in  mind,  we  now  present  algorithm  1,  the  basic  GJK 
algorithm.  The  algorithm  begins  by  setting  the  search  direction  to  be  the  x  direction 
and  creating  an  empty  simplex.  The  algorithm  next  enters  main  loop.  On  each  itera¬ 
tion,  GJK  uses  the  search  direction  to  find  a  new  point  on  the  convex  hull  of  Mdiff, 
and  adds  that  point  to  the  simplex.  If  the  origin  is  not  crossed  by  moving  to  this 
new  point,  then  GJK  reports  no  intersection.  Otherwise,  the  simplex  is  processed  to 
determine  if  the  origin  is  contained  inside.  If  so,  then  GJK  reports  that  the  2  objects 
intersect.  If  not,  then  GJK  gets  a  new  simplex  and  search  direction  to  use  in  the  next 
iteration.  We  discuss  how  to  process  the  simplex  in  Section  3.2. 

The  GJK  algorithm  terminates  when  either  of  the  terminating  conditions  discussed 
in  Section  3.1.3  are  met.  The  algorithm  also  terminates  after  20  iterations  have  been 
performed,  which  we  added  to  ensure  that  the  algorithm  would  always  terminate. 
The  value  of  20  is  somewhat  arbitrary.  In  our  experience,  GJK  generally  finds  an 
answer  within  a  handful  of  iterations  and  rarely  needs  more  than  10.  We  also  found 
that  if  GJK  got  to  20  iterations,  there  was  usually  a  problem  or  bug  in  the  code  that 
needed  to  be  solved.  For  now,  we  leave  the  iteration  cap  in  the  algorithm  to  ensure 
the  algorithm  always  terminates. 

The  top  level  algorithm  makes  use  of  a  few  additional  methods:  add(),  num_points(), 
processO,  contains_origin(),  and  get_next().  We  do  not  provide  any  pseudo-code 
for  these  methods  as  they  involve  book  keeping  details  that  detract  from  the  over¬ 
all  algorithm.  Nonetheless,  their  operations  should  be  clear.  add()  simply  stores  the 
point  in  whatever  data  structure  the  simplex  uses.  num_points()  returns  the  num¬ 
ber  of  points  in  the  simplex.  process()  calls  the  appropriate  process  methods  given 
below.  contains_origin()  returns  true  if  the  simplex  contains  the  origin.  get_next() 
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def  collision(,<2/,  ^): 
niters  =  0 
done  =  False 
res  =  False 
dir  =  [1,  0,  0] 
simp  =  SimplexO 
while  not  done  and  niters  <  20: 
niters  +=  1 

P  =  ■c/.support(dir)  — e^.support(— dir) 
simp.add(F) 

if  simp.num_points()  >  1  and  OP  ■  dir  <  0: 
res  =  False 
done  =  True 
simp.processO 
if  simp.contains_origin(): 
res  =  True 
done  =  True 

simp,  dir  =  simp.get_next() 

return  res 

Algorithm  1:  Basic  GJK  algorithm 


returns  the  simplex  and  direction  to  use  in  the  next  iteration  of  GJK. 

3.2  Simplex  Processing 

At  each  iteration  of  the  GJK  algorithm,  GJK  adds  a  point  to  the  current  simplex  and 
then  processes  the  simplex.  As  we  saw  in  the  examples  in  Section  3.1.2,  processing 
the  simplex  involves  determining  if  the  origin  is  inside  the  simplex.  If  so,  then  there 
is  no  further  action  required.  Otherwise,  GJK  examines  the  simplex  to  remove  any 
unnecessary  points  and  then  computes  a  new  direction  to  search  for  the  next  point 
to  add. 

We  determine  if  the  origin  is  inside  the  simplex  and  which  points  are  unnecessary 
through  the  same  general  process.  For  each  simplex,  we  examine  each  geometrical 
feature  (i.e.,  vertex  points,  line  segments,  and  triangle  faces)  one  by  one.  Some  fea¬ 
tures  can  immediately  be  eliminated  from  consideration  simply  due  to  the  way  the 
simplex  was  constructed.  For  each  remaining  feature,  we  compute  a  direction, 
d,  that  points  from  ^  and  points  away  from  the  interior  of  the  simplex.  Usually  d  is 
the  vector  normal  to  ^ ,  but  not  always.  We  next  compute  a  vector,  FO,  that  points 
from  to  the  origin.  If  d  has  a  positive  projection  along  FO,  i.e.,  FO  ■  d  >  0,  then 
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we  conclude  that  the  origin  lies  outside  the  simplex  and  that  ^  is  the  feature  that 
we  need.  Any  point  in  the  simplex  that  is  not  part  of  ^  is  no  longer  needed  and  can 
be  removed  the  simplex.  Thus,  the  simplex  for  the  next  iteration  simply  contains 
^ .  The  search  direction  for  the  next  iteration  is  d. 

If  no  such  feature  can  be  found,  then  the  origin  lies  within  the  simplex,  or  possibly 
on  one  of  the  features  themselves.  Regardless,  we  conclude  that  the  origin  is  con¬ 
tained  with  the  simplex,  which  means  that  we  have  discovered  that  the  2  objects 
intersect. 

Simplex  processing  can  be  performed  by  using  dot  and  cross  products  appropriately 
applied  to  each  feature,  which  makes  the  process  is  more  geometrical  in  nature.  By 
contrast,  the  original  GJK  paper^  determined  which  feature  was  closest  to  the  origin 
through  a  very  algebraic  procedure.  The  geometrical  approach  is  generally  simpler 
to  implement  and  more  intuitive  than  the  algebraic  approach.  However,  we  must  be 
extremely  careful  in  how  the  direction  d  is  computed  from  each  feature  ^ .  If  d  is 
computed  incorrectly,  then  GJK  may  not  be  able  to  determine  if  the  origin  is  located 
within  a  simplex.  A  number  of  the  bugs  we  discovered  during  the  development 
process  were  due  to  slight  mistakes  in  computing  the  directions. 

We  also  remark  that  the  search  direction  computed  for  the  next  iteration  of  GJK  is 
always  from  a  feature  and  not  from  a  point.  (The  only  exception  being  the  second 
iteration  of  GJK  because  the  first  iteration  creates  a  0-simplex.)  This  can  have  some 
surprising  results,  as  we  saw  in  Section  3.1.2  when  examining  the  left  side  of  Fig.  4. 
After  we  added  point  E  to  the  simplex,  we  next  added  point  G  and  not  F ,  despite 
the  fact  that  F  would  appear  to  be  closer  to  the  origin  and  E.  The  reason  for  this  is 
that  the  direction  computed  was  from  the  line  segment  feature  AE  and  not  the  point 
E  itself.  We  discuss  this  more  when  we  handle  the  1 -simplex  below. 

When  discussing  each  simplex  type,  we  label  the  points  in  the  simplex  as  A,  B,  C, 
and  D.  We  adopt  the  convention  that  A  is  always  the  point  that  was  most  recently 
added  to  the  simplex.  This  convention  is  critical,  as  it  allows  us  to  skip  certain 
checks  that  we  would  otherwise  need  to  perform.  The  ordering  of  the  other  points 
is  less  critical,  though  we  generally  make  B  the  second  most  recently  added  point, 
C  the  third  most  recently  added  point,  and  D  the  first  point  added.  As  we  discuss 
below,  there  are  2  instances  where  this  ordering  is  modified:  in  the  2-simplex  pro¬ 
cessing  and  in  the  3-simplex.  Regardless,  A  is  always  the  most  recently  added  point. 
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Finally,  for  each  simplex  type,  we  present  pseudo-code  for  processing  that  simplex. 
The  pseudo-code  uses  the  notation  given  in  Section  2,  e.g.,  AB  refers  to  the  vector 
from  point  A  to  point  B,  AB  x  AC  is  the  vector  cross  product,  and  so  on.  The  syntax 
of  the  pseudo-code  is  an  approximate  object-oriented  syntax  where  “self”  refers  to 
the  instance  variable  (C-i-i-  programmers  would  use  “this”).  We  present  the  pseudo¬ 
code  for  processing  each  simplex  as  methods  of  the  simplex  class. 

3.2.1  Processing  a  0-simplex 

A  0-simplex  consists  of  a  single  point.  A,  that  lies  on  the  convex  hull  of  Mdiff.  Since 
this  vertex  point  is  the  only  feature  in  the  simplex,  there  are  no  further  checks  to  be 
performed.  The  simplex  for  the  next  iteration  contains  this  point.  The  only  direction 
we  can  move  is  back  from  A  toward  the  origin.  Since  A  is  the  vector  from  the  origin 
to  A,  —A  must  be  the  direction  we  need  to  move  to  come  back  to  the  origin. 

A  very  rare  case,  although  one  that  may  still  happen,  is  that  point  A  may  actually 
be  the  origin.  In  order  for  the  origin  to  lie  on  the  exterior  of  Mdiff,  the  2  objects 
must  intersect  at  exactly  1  point.  Whether  this  contact  constitutes  an  intersection 
depends  on  the  application.  Some  applications  regard  such  contact  as  a  collision 
while  others  do  not.  We  regard  such  contact  as  a  collision,  and,  thus,  need  to  check 
for  the  case  where  A  is  the  origin.  For  applications  that  do  not  consider  this  case 
to  be  an  intersection,  the  algorithm  for  processing  the  0-simplex  would  need  to  be 
adjusted,  as  would  some  portions  of  the  steps  for  processing  the  other  simplices. 
With  this  in  mind,  we  present  algorithm  2  to  handle  the  processing  of  a  0-simplex. 


def  process_Osimplex(self): 
if  A  is  the  origin: 

mark  self  as  containing  the  origin 

else: 

self.next_direction  =  —A 
self.next_simplex  =  self 

Algorithm  2:  Processing  a  0-simplex 


3.2.2  Processing  a  1-simplex 

A  1-simplex  consists  of  the  line  segment  between  the  2  points.  A,  and  B,  where 
A  was  the  most  recent  point  added.  There  are  3  features,  the  point  A,  the  point  B, 
and  the  line  segment  AB.  These  features  divide  the  space  around  the  simplex  into 
4  different  regions,  as  shown  in  Fig.  6.  Region  Ri  lies  beyond  B.  Region  R2  lies 
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Fig.  6  Sketch  of  a  1-simplex,  which  is  a  line  segment.  The  origin  could  he  in  any  of  the  4 
indicated  regions. 


beyond  A.  Regions  Rt,  and  R4  lie  on  either  side  of  the  line  segment.  The  origin 
could  lie  in  any  of  these  regions  or  on  the  line  segment  itself. 

We  need  to  check  each  of  the  regions  for  the  origin.  We  can  immediately  eliminate 
Ri  and  R2  from  consideration.  At  the  beginning  of  this  iteration,  GJK  had  a  0- 
simplex  containing  the  point  B.  GJK  moved  away  from  B  toward  the  origin  and 
found  the  point  A.  Thus,  the  origin  cannot  be  in  region  Ri  since  GJK  moved  away 
from  that  region  already.  The  origin  also  cannot  be  in  R2.  If  the  origin  was  in  R2,  the 
GJK  would  have  terminated  since  going  from  BtoA  did  not  cross  over  the  origin. 

That  means  the  origin  is  located  in  either  R3  or  R4.  We  handle  both  simultaneously. 
AB  is  the  vector  from  A  to  5  while  AO  is  the  vector  from  A  to  the  origin.  Those  2 
vectors  define  a  plane  whose  normal  is  given  by  AB  x  AO.  If  we  cross  the  normal 
into  AB  again,  we  obtain  a  vector  that  points  from  the  line  segment  toward  the 
origin.  Thus,  the  direction  we  want  is  given  by  AB  x  AO  x  AB.  Note  that  it  does 
not  matter  if  the  origin  lies  in  R3  or  R4  since  both  cases  result  in  the  same  expression 
for  the  direction. 

Here  we  could  ask  why  choose  AB  x  AO  x  AB  as  the  new  direction  to  check,  as 
opposed  to  simply  choosing  AO  since  AO  points  directly  to  the  origin.  The  reason 
is  that  AB  x  AO  x  AB  is  perpendicular  to  the  line  segment,  but  AO  may  not  be.  In 
fact,  AO  could  be  located  very  close  to  AB.  If  we  move  in  that  direction,  we  could 
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obtain  B  as  our  point  to  use  for  the  next  iteration.  This  would  cause  a  ping-pong  like 
effect  where  GJK  continually  bounces  between  selecting  B,  then  A,  then  B,  then  A, 
and  so  on.  Since  we  determined  that  the  origin  lies  closest  to  the  line  segment  and 
not  either  of  the  vertex  points,  we  require  a  direction  from  the  line  segment  towards 
the  origin.  The  only  useful  direction  would  be  AB  x  AO  x  AB. 

We  also  need  to  consider  2  boundary  cases:  A  could  be  the  origin  itself  or  the  origin 
could  lie  on  the  line  segment  A5. 

•  If  A  is  the  origin,  then  AO  is  zero,  meaning  that  AB  x  AO  x  AB  =  0. 

•  If  the  origin  lies  on  the  segment  AB,  then  AB  x  AO  =  0  since  the  2  vectors 
would  be  parallel. 

Thus,  both  of  these  cases  require  checking  if  AB  x  AO  x  AB  is  zero.  Some  imple¬ 
mentations  do  not  check  for  these  cases. Not  handling  these  cases  forces  GJK 
to  use  a  direction  of  zero  to  search  for  the  next  point  to  add.  This  could  cause  prob¬ 
lems  when  used  in  the  support  functions  since  support  functions  are  to  return  a  point 
maximally  far  in  the  given  direction.  The  point  maximally  far  in  the  zero  direction 
is  an  undefined  concept,  so  the  support  functions  may  act  in  unknown  ways  when 
given  a  zero  vector.  An  example  of  this  is  the  support  function  for  a  sphere  given 
in  Eq.  6.  If  given  a  zero  direction,  that  support  function  would  return  the  center  of 
the  sphere  and  not  a  point  on  the  edge  of  the  sphere.  This  would  cause  GJK  to  con¬ 
sider  points  that  are  not  on  the  exterior  of  Mdiff,  which  violates  our  assumptions. 
We  could  fix  this  behavior  by  demanding  the  support  functions  check  for  a  zero 
direction,  but  we  feel  it  is  better  to  simply  handle  the  case  in  our  implementation  of 
GJK. 

def  process_lsimplex(self): 

dir  =  AB  x  AO  x  AB 
if  dir  is  zero: 

mark  self  as  containing  the  origin 

else: 

self.next_direction  =  dir 
self.next_simplex  =  self 

Algorithm  3:  Processing  a  1- simplex 

Algorithm  3  shows  our  implementation.  Since  the  line  segment  AB  is  the  feature 
closest  to  the  origin,  we  require  both  points  in  the  new  simplex. 
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3.2.3  Processing  a  2-simplex 


Fig.  7  Sketch  of  a  2-simplex  and  the  regions  created  around  it. 


A  2-simplex  is  considerably  more  complex  to  handle  than  either  the  1 -simplex  or 
0-simplex.  The  2-simplex  forms  a  triangle  in  space,  which  contains  7  features:  3 
vertex  points,  3  line  segments,  and  the  triangle  face  itself.  Each  of  these  features 
corresponds  to  a  particular  region  in  space,  as  shown  in  Fig.  7.  We  label  the  vertices 
of  the  triangle.  A,  B,  and  C,  again  with  the  understanding  that  A  was  the  last  point 
added  to  the  simplex.  The  7  regions  are  Ri,  which  is  the  interior  of  the  triangle;  R2, 
R4  and  R(i,  which  are  the  regions  beyond  the  vertices;  and  R^,  R5  and  R^,  which  lay 
beyond  the  line  segments  between  each  parr  of  vertices. 

We  now  need  to  determine  which  region  the  origin  lies  in.  Of  the  7  possible  regions, 
we  immediately  eliminate  3  from  consideration.  In  the  previous  iteration,  GJK  had 
a  1 -simplex  consisting  of  the  points  B  and  C.  GJK  sought  a  new  point  by  searching 
from  BC  toward  the  origin,  resulting  in  point  A.  Thus,  in  searching  for  a  point  toward 
the  origin,  GJK  moved  away  from  R4,  Rs  and  Rf,,  so  the  origin  cannot  be  in  those 
regions. 

We  can  also  eliminate  region  R2  from  consideration.  If  the  origin  were  in  R2,  then 
moving  from  BC  to  point  A  would  not  have  crossed  the  origin,  which  means  the  top 
level  GJK  algorithm  would  have  terminated  already.  The  origin  therefore  must  be 
located  in  regions  R^,  Rj,  or  Ri. 

Since  R3  and  R^  are  regions  that  lie  beyond  line  segments,  we  would  expect  their 
handling  to  be  similar  to  that  of  a  1 -simplex.  Once  again,  we  define  a  normal  to  the 
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line  segment  and  check  the  region  for  the  origin.  However,  there  are  important  dif¬ 
ferences  from  the  1 -simplex  case  since  A5  and  AC  come  from  a  triangle.  We  cannot 
simply  use  the  same  procedure  from  1 -simplex  case  because  we  must  differentiate 
between  the  cases  of  the  origin  lying  in  R^,  or  or  Ri. 

Our  task  is  to  define  a  normal  from  AB  that  points  into  Rt,  and  a  normal  from  AC 
that  points  into  R-j,  and  then  check  each  region.  In  3  dimensions,  there  are  an  infinite 
number  of  vectors  normal  to  either  line  segment.  However,  only  the  normals  in  the 
plane  of  the  triangle  are  useful.  Each  segment  has  2  such  (unit)  normals,  pointing 
outward  from,  and  inward  to  7?i.  We  determine  these  normals  by  first  computing 
ABC„,  which  is  the  normal  to  the  triangle, 

ABC„=ABxAC  (10) 

The  normals  to  each  of  the  line  segments  are  then  defined  as 

AB„  =  AB  X  ABC„  AC„  =  ABC„  x  AC  (1 1) 

Carefully  note  the  order  in  which  the  vector  cross  products  are  computed.  We  want 
the  normals  to  the  line  segments  to  point  away  from  the  interior  of  the  triangle.  The 
normals  defined  in  Eq.  1 1  obey  this  convention.  If  the  order  of  the  vectors  are  flipped 
in  the  cross  products,  then  the  resulting  normals  will  point  toward  the  interior  of  the 
triangle,  which  will  cause  GJK  to  compute  incorrect  results.  Our  experience  is  that 
computing  incorrect  normals  is  a  common  way  of  introducing  bugs  into  the  code, 
which  can  be  very  difficult  to  find. 

With  the  normals  correctly  defined,  we  determine  if  the  origin  lies  in  Rt,  by  checking 
the  condition 

AB„-AO>0.  (12) 

If  Eq.  12  is  true,  then  the  origin  lies  somewhere  in  R^,  so  we  no  longer  need  point  C 
in  our  simplex.  We  now  need  a  direction  to  search  for  the  next  point.  We  compute 
that  direction  as 

ABxAOxAB,  (13) 

which  is  analogous  to  the  direction  used  for  the  1 -simplex  case.  Note  that  the  di¬ 
rection  we  need  to  move  is  not  given  by  AB„.  In  2  dimensions,  we  could  use  AB„ 
as  our  new  direction,  but  in  3  dimensions  the  origin  may  lie  above  (or  below)  the 
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plane  of  the  triangle.  The  direetion  in  Eq.  13  takes  this  into  aeeount,  whereas  AB„ 
does  not.  The  eheek  for  is  identieal,  exeept  that  C  is  used  instead  of  B. 

For  eompleteness  in  handling  and  R^,  we  note  that  we  found  a  referenee^  that 
advoeated  eheeking  AB  ■  AO  >  0  and  BA  ■  BO  >  0  in  addition  to  Eq.  12  for  R^  (and 
similar  for  R^).  In  general,  we  do  not  feel  these  extra  eonditions  are  neeessary  as 
they  are  almost  always  satisfied.  In  order  for  the  2-simplex  to  be  proeessed,  we  must 
have  erossed  the  origin  when  moving  to  point  A,  as  per  the  termination  eondition 
in  GJK  algorithm  (algorithm  1).  That  guarantees  that  AB  ■  AO  >  0  if  A  is  the  point 
exaetly  on  the  direetion  we  moved  from  B.  We  note  that  the  support  funetions  do 
allow  for  some  “wiggle”  room  in  that  they  do  not  have  to  report  a  point  exaetly 
in  the  direetion.  However,  the  support  funetions  still  return  a  point  approximately 
in  the  direetion,  whieh  is  usually  enough  to  avoid  these  extra  eheeks.  Nonetheless, 
there  may  be  some  peeuliar  eases  that  require  these  extra  eheeks,  though  we  never 
eneountered  any  in  our  experienee  with  GJK.  Other  referenees^"^  omit  the  extra 
eheeks  entirely. 

If  the  origin  is  not  loeated  within  R^  or  Rj,  then  it  must  be  loeated  within  Ri,  whieh 
is  the  interior  of  the  triangle.  In  2  dimensions,  this  would  be  enough  to  report  that 
the  2  objeets  interseet.  However,  in  3  dimensions,  the  origin  may  lie  in,  above,  or 
below  the  plane  of  the  triangle,  eaeh  of  whieh  needs  to  be  handled  separately. 

•  If  ABC„  ■  AO  =  0,  then  the  origin  is  in  the  plane  of  the  triangle  and  loeated  in 
the  interior  of  the  triangle.  This  means  we  have  found  the  origin  within  Mdiff 
and  ean  report  an  interseetion. 

•  If  ABC„  ■  AO  >  0,  then  the  origin  lies  above  the  plane  of  the  triangle.  In 
this  ease,  the  simplex  stays  the  same  and  we  seareh  for  the  next  point  in  the 
direetion  of  ABC„. 

•  If  ABC„  ■  AO  <  0,  then  the  origin  lies  below  the  plane  of  the  triangle.  We  thus 
need  to  move  in  the  direetion  of  — ABC„  to  find  our  next  point.  Unlike  the 
earlier  ease  where  the  origin  lay  above  the  plane,  we  need  to  adjust  the  sim¬ 
plex  slightly.  As  diseussed  in  the  next  seetion,  the  proeessing  of  a  3-simplex 
requires  that  AB  x  AC  point  in  the  direetion  of  the  next  point  added.  This 
is  not  true  if  we  move  below  the  plane  of  the  triangle.  To  make  it  true,  we 
interehange  A  and  B  in  our  simplex  ordering. 
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We  also  need  to  eonsider  the  possibility  that  A  might  be  the  origin,  or  that  the  origin 
might  lie  on  1  of  the  line  segments.  As  it  turns  out,  none  of  these  eases  requires  any 
speeial  handling. 

•  If  A  is  the  origin,  then  AO  will  be  zero,  whieh  means  all  the  dot  produets 
involving  AO  will  also  be  zero.  By  the  above  diseussion,  the  origin  would  lie 
in  7?i,  and,  sinee  ABC„  ■  AO  =  0,  the  origin  would  lie  within  the  plane  of  the 
triangle.  We  thus  report  an  interseetion. 

•  If  the  origin  lies  on  any  of  the  edges,  then  the  dot  produet  of  AO  with  the 
normal  to  that  edge  will  be  zero.  The  dot  produet  with  the  other  edge  normal 
will  be  negative,  so  the  origin  again  lies  in  Ri.  Sinee  the  origin  is  in  the  plane 
of  the  triangle,  we  report  an  interseetion. 

With  these  diseussions  in  mind,  the  algorithm  for  handling  the  2-simplex  ease  is 
shown  in  algorithm  4. 

3.2.4  Processing  a  3-simplex 

A  3-simplex  is  a  tetrahedron  eonsisting  of  4  triangle  faees.  As  sueh,  we  would 
expeet  that  proeessing  a  3-simplex  is  similar  to  proeessing  a  2-simplex.  There  are 
still  a  few  slight  differenees,  however. 

Figure  8  shows  a  3-simplex,  whieh  is  a  useful  diagram  for  the  rest  of  this  seetion. 
GJK  started  the  previous  iteration  with  a  2-simplex  that  eonsisted  of  the  points  B, 
C  and  D.  During  that  iteration,  point  A  was  added  to  the  simplex  to  ereate  the  3- 
simplex  that  we  now  need  to  examine.  Importantly,  we  make  the  assumption  that 
the  point  A  lies  in  the  direetion  of  BC  x  BD.  If  this  assumption  is  violated,  then  our 
proeessing  of  the  3-simplex  will  run  into  diffieulties. 

To  ensure  this  assumption  is  met,  we  need  to  reexamine  how  the  3-simplex  was 
ereated  from  a  2-simplex.  During  the  proeessing  of  the  2-simplex,  GJK  ehose  the 
next  seareh  direetion  to  be  above  or  below  the  plane  of  the  triangle  to  find  point  A.  If 
the  direetion  was  above  the  plane,  then  the  assumption  of  BC  x  BD  pointing  toward 
A  is  automatieally  satisfied.  If  the  direetion  was  below  the  plane  of  the  triangle,  then 
A  points  in  the  direetion  of  — BC  x  BD,  whieh  violates  our  assumption.  However, 
if  we  exehange  the  order  of  any  2  of  the  points,  say  B  and  C,  then  the  direetion 
ehanges.  This  is  the  reason  we  ehanged  the  order  of  the  points  in  the  2-simplex 
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def  process_2simplex(self): 

ABC„  =  AB  X  AC 
AC„  =  ABC„  X  AC 
AB„  =  AB  X  ABC„ 
if  AB„  AO>0: 

#  Region  3 

self.next_direction  =  AB  x  AO  x  AB 
self.next_simplex  =  Simplex([B,  A]) 
elif  AC„  AO>0: 

#  Region  7 

self.next_direction  =  AC  x  AO  x  AC 
self.next_simplex  =  Simplex([C,  A]) 

else: 

#  Region  1;  Origin  could  be  below,  or 

#  above,  or  in  the  plane  of  the  triangle 

v  =  ABC„  AO 

if  V  is  zero: 

#  In  the  plane  of  the  triangle 

mark  self  as  containing  the  origin 

elifv>0: 

#  Above  the  plane  of  the  triangle 

self.next_direction  =  ABC„ 
self.next_simplex  =  self 

else: 

#  Below  the  plane  of  the  triangle 

#  Note:  we  change  the  order  of  the  simplex.  See 
the 

#  handling  of  a  3-simplex  for  why. 
self.next_direction  =  — ABC„ 
self.next_simplex  =  Simplex([C,  A,  B]) 

Algorithm  4:  Processing  a  2- simplex 
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Fig.  8  Sketch  showing  a  3-simplex  (a  tetrahedron).  The  point  A  lies  above  the  plane  deter¬ 
mined  hy  the  triangle  ABCD.  The  point  B  lies  behind  the  AACD  triangle  face. 


when  the  search  was  chosen  to  point  below  the  triangle.  By  switching  2  of  the 
points,  we  inverted  the  direction  of  the  normal  and  ensured  this  assumption  is  valid. 

With  our  tetrahedron  properly  oriented,  we  now  need  to  determine  which  region  the 
origin  could  lie  in.  We  have  15  features  to  examine:  4  vertices,  6  line  segments,  4 
triangle  faces,  and  the  tetrahedron  itself.  As  we  have  done  previously,  we  imme¬ 
diately  eliminate  some  features  from  consideration.  The  previous  iteration  used  a 
2-simplex  that  had  the  points  B,  C,  and  D,  and  moved  toward  the  origin  to  find  the 
point  A.  Thus,  we  can  eliminate  the  regions  related  to  the  points  B,  C,  and  D,  the 
line  segments  BC,  CD  and  DB,  and  the  triangle  ABCD.  Further,  we  can  eliminate 
the  region  beyond  A.  If  the  origin  was  there,  then  GJK  would  have  terminated  al¬ 
ready  since  moving  from  ABCD  to  A  did  not  cross  the  origin.  This  leaves  us  with 
7  features:  the  tetrahedron,  3  triangle  faces,  and  3  line  segments. 

We  process  those  features  by  first  determining  if  the  origin  lies  above  any  of  the 
triangle  faces,  AABC,  AACD,  and  AADB.  If  so,  then  we  conclude  that  the  origin  is 
outside  the  tetrahedron,  but  may  possibly  still  lie  within  Mdiff.  Thus,  a  new  simplex 
must  be  constructed  and  a  new  search  direction  chosen.  The  new  simplex  is  either 
the  triangle  face  itself,  or  an  edge,  whichever  feature  is  closer  to  the  origin.  If  the 
origin  does  not  lie  outside  of  any  of  the  faces,  then  the  origin  is  contained  within 
the  tetrahedron  and  we  can  report  that  the  objects  intersect. 
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We  determine  if  the  origin  lies  above  a  triangle  faee  by  eomputing  the  dot  produet 
of  AO  with  the  normal  to  the  triangle  faee.  If  the  dot  produet  is  positive,  then  the 
origin  lies  elosest  to  that  faee.  For  this  proeedure  to  work,  we  need  to  define  the 
normals  to  eaeh  faee  so  that  the  normals  point  outward  from  the  tetrahedron.  We 
define  the  normals  as 

ABC„  =  ABxAC  ACD„=ACxAD  ADB„  =  AD  x  AB.  (14) 


Note  earefully  how  the  normals  are  defined  in  Eq.  14.  The  order  of  the  veetors  in  the 
eross  produets  were  ehosen  to  ensure  the  normals  pointed  outward.  Choosing  any 
other  order  eauses  at  least  1  of  the  normals  to  point  inward.  Further,  the  definitions 
of  the  normals  are  eonsistent  with  the  tetrahedron  laid  out  in  Fig.  8.  This  is  why  we 
made  the  assumption  of  A  lying  in  the  direetion  of  BC  x  BD.  If  A  does  not  lie  in 
that  direetion,  then  normals  would  have  to  be  defined  differently. 

We  now  have  to  determine  how  to  handle  eaeh  triangle  faee.  The  proeedure  is  iden- 
tieal  for  eaeh  faee,  so  we  only  need  to  diseuss  the  proeess  for  a  single  faee.  Suppose 
that  AO  ■  ABC„  >  0,  then  the  origin  lies  elosest  to  the  AABC  triangle.  Right  away, 
we  know  that  point  D  is  no  longer  needed  and  ean  be  removed  from  the  simplex.  We 
now  need  to  determine  if  the  origin  lies  eloser  to  an  edge  or  to  the  faee  itself.  This 
is  almost  identieal  to  the  2-simplex  handling  diseussed  earlier.  The  only  differenee 
is  that  we  do  not  need  to  eheek  if  the  point  lies  below  the  triangle  sinee  we  already 
know  that  it  lies  above  it.  That  leaves  us  with  3  possibilities:  the  origin  lies  in  the 
direetion  from  the  AC  edge  or  the  AB  edge,  or  above  the  AABC  faee  itself. 

As  with  the  2-simplex  ease,  we  distinguish  between  these  eases  by  eomputing  the 
relevant  normals.  For  the  edges,  we  define  the  normals  as  AC^  =  ABC^  x  AC  and 
AB„  =  AB  X  ABC„.  The  triangle  faee  has  normal  ABC„,  whieh  we  already  eom- 
puted.  We  eheek  for,  and  handle,  eaeh  of  the  3  eases  by  eheeking  the  dot  produet  of 
the  normal  with  AO. 

•  If  AO  ■  AC„  >  0,  then  the  origin  lies  eloser  to  the  AC  edge.  We  no  longer  need 
points  B  and  D,  so  we  ereate  a  1 -simplex  eonsisting  of  A  and  C.  The  new 
direetion  to  move  is  AC  x  AO  x  AC. 

•  If  AO  ■  AB„  >  0,  then  the  origin  lies  eloser  to  the  AB  edge.  This  is  identieal 
to  the  first  ease,  exeept  that  B  and  C  interehange. 
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•  If  neither  of  those  cases  is  true,  then  the  origin  lies  directly  above  the  triangle 
face.  We  could  check  AO  ■  ABC„  to  confirm  this,  but  we  already  computed 
this  product.  In  this  case,  we  create  a  new  2-simplex  consisting  of  A,  B,  and 
C.  The  new  direction  to  search  is  given  by  ABC„. 

The  processing  of  the  triangles  AACD  and  AADB  is  almost  identical  to  the  AABC 
case  except  for  changing  which  points  are  involved.  There  are,  however,  2  caveats 
that  need  careful  attention.  The  first  caveat  is  that  the  edge  normals  are  computed 
differently  when  using  different  triangle  normals.  Thus,  AC„  is  different  when  ex¬ 
amining  triangle  AABC  than  it  is  when  examining  triangle  AACD. 

The  second  caveat  is  that  when  creating  a  2-simplex,  we  need  to  ensure  that  the 
points  are  ordered  correctly  in  the  simplex.  If  we  decide  to  create  a  2-simplex, 
then,  in  the  next  iteration,  we  add  another  point  creating  a  3-simplex.  Consistent 
with  our  assumption  above,  we  require  that  new  point  to  lie  in  the  direction  of  the 
cross  products  of  the  edges.  Thus,  we  need  to  order  the  points  to  ensure  that  this 
assumption  is  upheld.  This  only  affects  the  case  for  the  AADB  triangle. 

Finally,  the  origin  could  either  be  point  A  itself,  or  the  origin  could  lie  in  1  of  the 
triangle  faces,  or  the  origin  could  lie  in  the  plane  of  1  of  the  triangles,  but  not  in  the 
face  itself.  These  boundary  cases  do  not  require  any  special  processing. 

•  If  A  is  the  origin,  then  AO  will  be  zero  and  all  the  dot  products  we  check 
will  be  zero.  Thus,  we  consider  that  the  origin  lies  inside  the  tetrahedron  and 
report  an  intersection. 

•  If  the  origin  lies  within  a  triangle  face,  then  the  dot  product  of  AO  with  that 
face’s  normal  will  be  zero  and  the  dot  product  with  the  other  normals  will 
be  negative.  Again,  this  causes  us  to  consider  the  origin  as  being  inside  the 
tetrahedron. 

•  Lastly,  if  the  origin  lies  in  the  plane  of  a  triangle,  but  not  in  the  triangle  face 
itself,  then  the  dot  product  of  AO  with  the  plane  normal  will  still  be  zero.  In 
this  case,  the  origin  will  lies  above  1  of  the  other  triangle  faces  causing  the 
dot  product  of  AO  with  that  normal  to  be  positive,  and  so  this  case  is  handled 
by  creating  a  new  2-simplex  and  moving  on. 

With  that  in  mind,  the  algorithm  for  handling  the  3-simplex  is  shown  in  algorithm  5. 
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def  process_3simplex(self): 

#  The  3-simplex  consists  of  the  points  A,  B,  C,  D. 

#  We  assume  that  A  lies  in  the  direction  of  BC  x  BD 

ABC„  =  AB  X  AC 
ACD„  =  AC  X  AD 
ADB„  =  AD  X  AB 
if  ABC„  AO>0: 

AC„  =  ABC„  X  AC 
AB„  =  AB  X  ABC„ 
if  AC„  AO>0: 

self.next_simplex  =  Simplex([C,  A]) 
self.next_direction  =  AC  x  AO  x  AC 
elif  AB„  ■  AO  >  0: 

self.next_simplex  =  Simplex([B,  A]) 
self.next_direction  =  AB  x  AO  x  AB 
else: 

self.next_simplex  =  Simplex([C,  B,  A]) 
self.next_direction  =  ABC„ 
elif  ACD„  AO>0: 

AD„  =  ACD„  X  AD 
AC„  =  AC  X  ACD„ 
if  AD„  AO>0: 

self.next_simplex  =  Simplex([D,  A]) 
self.next_direction  =  AD  x  AO  x  AD 
elifAC„  AO>0: 

self.next_simplex  =  Simplex([C,  A]) 
self.next_direction  =  AC  x  AO  x  AC 
else: 

self.next_simplex  =  Simplex([D,  C,  A]) 
self.next_direction  =  ACD„ 
elif  ADB„  AO>0: 

AB„  =  ADB„  X  AB 
AD„  =  AD  X  ADB„ 
if  AB„  AO>0: 

self.next_simplex  =  Simplex([B,  A]) 
self.next_direction  =  AB  x  AO  x  AB 
elif  AD„  ■  AO  >  0: 

self.next_simplex  =  Simplex([D,  A]) 
self.next_direction  =  AD  x  AO  x  AD 
else: 

self.next_simplex  =  Simplex([B,  D,  A]) 
self.next_direction  =  ADB„ 

else: 

mark  self  as  containing  the  origin 

Algorithm  5:  Algorithm  for  processing  a  3-simplex 
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3.3  GJK  Efficiency;  Qualitative 


GJK  tends  to  be  extremely  efficient,  both  spatially  and  temporarily,  when  com¬ 
puting  answers.  There  are  no  large  data  structures  to  process,  nor  any  particularly 
computationally  intensive  steps  in  the  algorithm.  These  factors  make  GJK  quite  ef¬ 
ficient.  Establishing  absolute  bounds  on  the  spatial  and  temporal  requirements  is 
difficult  because  GJK  relies  upon  external  data  structures  (the  objects  being  tested) 
and  external  support  functions.  These  influences  could  cause  GJK  to  vary  wildly 
in  how  much  space  and  time  are  required  to  test  the  intersection  of  2  objects.  Pro¬ 
vided  that  these  external  influences  are  not  too  burdensome,  GJK  performs  very 
efficiently. 

Spatially,  GJK  only  requires  space  for  the  simplex  (ignoring  the  space  for  the  ob¬ 
jects).  The  simplex  only  needs,  at  most,  four  3-D  points;  a  new  direction  to  search 
for  points;  a  few  state  variables,  such  as  a  flag  describing  whether  the  simplex  con¬ 
tains  the  origin;  and  space  for  a  few  temporary  vectors  used  during  the  processing  of 
each  simplex  type.  None  of  these  items  are  large,  making  the  spatial  requirements 
almost  negligible. 

Temporally,  GJK  requires  a  series  of  iterations  to  determine  the  result.  The  single 
most  expensive  operation  in  any  of  the  iterations  is  computing  vector  cross  prod¬ 
ucts,  which  requires  6  multiplications  and  3  subtractions.  Processing  a  3-simplex 
requires,  at  most,  7  cross  products  to  be  computed,  which  does  not  require  much 
time  on  modern  processors. 

The  number  of  iterations  GJK  requires  is  difficult  to  predict,  as  it  depends  on  the 
objects  involved  and  the  initial  direction  chosen.  GJK  can  usually  figure  out  an 
answer  within  a  handful  of  iterations,  usually  around  5-10.  If  the  objects  intersect, 
GJK  generally  requires  at  least  4  iterations  in  order  to  create  the  3-simplex  and 
determine  if  the  origin  is  contained  within.  GJK  may  require  a  few  more  iterations 
to  construct  another  3-simplex,  but  usually  that  is  the  most  required.  If  the  objects 
do  not  intersect,  then  GJK  is  usually  able  to  determine  that  within  2-3  iterations 
since  GJK  quickly  discovers  that  the  origin  was  not  crossed  when  selecting  a  new 
point.  Thus,  we  generally  expect  no  more  than  about  10  iterations  to  determine  an 
answer. 

Algorithm  1  limited  the  number  of  iterations  to  20  to  ensure  that  the  algorithm 
always  terminates.  Twenty  was  arbitrarily  chosen  based  on  simply  doubling  the 
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number  of  expected  iterations.  In  our  experience,  we  found  that  the  only  time  GJK 
ever  needed  20  iterations  was  when  a  bug  was  present  in  the  code.  After  correcting 
any  bugs  that  occurred,  we  never  found  GJK  to  require  20  iterations.  We  left  the 
iteration  bound  in  the  algorithm  just  in  case  an  unknown  bug  was  encountered, 
although  it  may  be  possible  for  some  exotic  objects  to  require  that  many  iterations 
to  determine  collisions. 

Taken  all  together,  we  find  that  GJK  is  extremely  efficient  in  time  and  space  require¬ 
ments.  This  makes  GJK  applicable  in  a  wide  variety  of  applications,  especially  the 
system  effectiveness  modeling  that  spurred  our  initial  interest  in  GJK. 

3.4  GJK  Summary 

GJK  is  an  efficient  algorithm  for  detecting  the  collision  between  2  objects,  and 
SS.  GJK  reduces  the  intersection  problem  to  the  problem  of  determining  if  the  ori¬ 
gin  lies  within  Mdiff(.2/,  GJK  constructs  a  simplex  consisting  of  points  chosen 
from  the  convex  hull  of  Mdiff.  Each  simplex  divides  the  space  bounded  by  Mdiff 
into  several  regions  based  on  the  geometrical  features  (vertices,  line  segments,  tri¬ 
angle  faces)  of  the  simplex.  Each  region  is  checked  to  see  if  the  origin  may  lie  in 
that  region.  If  so,  then  GJK  creates  a  new  simplex  consisting  of  the  points  involved 
in  the  given  feature,  and  then  chooses  a  new  direction  to  search  for  a  new  point  from 
the  convex  hull  of  Mdiff.  The  direction  chosen  is  usually  perpendicular  to  feature 
itself.  If  no  region  contains  the  origin,  then  the  origin  lies  within  the  simplex  itself 
and  GJK  reports  an  intersection  since  Mdiff  would  therefore  contain  the  origin.  Of 
course,  the  origin  may  lie  on  a  feature  itself,  which  we  treat  as  special  cases.  If  GJK 
selects  a  new  point,  but  does  not  cross  over  the  origin  to  reach  that  point,  then  GJK 
stops  since  the  origin  is  clearly  not  contained  within  Mdiff. 

Algorithm  1  shows  the  basic  GJK  algorithm,  with  algorithms  2-5  showing  how 
to  process  each  type  of  simplex.  The  algorithms  show  that  GJK  is  fairly  simple  and 
straightforward  to  implement.  However,  as  the  discussion  about  those  algorithms 
indicates,  numerous  assumptions  and  conditions  go  into  the  algorithms  that  need 
careful  attention.  These  assumptions  and  conditions  are  not  often  stated,  which  can 
result  in  hard  to  find  bugs  when  implementing  the  algorithms.  Those  assumptions 
and  conditions  are  summarized  below: 

•  The  support  function  always  returns  a  point  farthest  in  the  given  direction. 

This  is  undefined  when  a  zero  direction  is  used,  so  we  need  to  ensure  that  the 
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support  functions  are  never  given  a  zero  direetion. 

•  The  points  added  to  the  simplex  always  eome  from  the  eonvex  hull  of  Mdiff. 
This  eondition  is  satisfied  provided  the  support  funetion  operates  eorreetly. 

•  In  the  algorithms,  the  point  A  is  always  the  last  point  added  to  the  simplex. 
The  ordering  of  the  other  points  in  the  simplex  is  less  eritieal,  but  we  generally 
maintain  the  order  in  whieh  they  were  added. 

•  When  we  proeess  the  2-simplex  and  3-simplex,  we  need  to  eompute  normals 
that  always  point  away  from  the  simplex.  The  normals  should  never  point 
toward  the  interior  of  the  simplex. 

•  For  a  3-simplex,  the  point  A  always  lies  in  the  direetion  of  BC  x  BD.  Obeying 
this  assumption  may  require  slightly  reordering  the  points  in  the  simplex. 

4.  Verification  of  Implementation 

After  implementing  the  GJK  algorithm,  we  sought  to  verify  our  implementation  by 
eomparing  results  against  an  algorithm  that  eomputes  interseetions  among  trian¬ 
gles,^®  whieh  we  label  TriTri.  TriTri  shares  nothing  in  eommon  with  GJK.  We  did 
not  ereate  our  own  implementation  TriTri,  ehoosing  instead  to  use  the  implementa¬ 
tion  ereated  by  TriTri ’s  authors.^® 

To  make  the  eomparison,  we  randomly  generated  1  million  pairs  of  triangles  in 
3  dimensions.  For  eaeh  pair,  we  used  both  our  GJK  implementation  and  TriTri  to 
determine  if  the  2  triangles  interseeted.  We  eounted  the  number  of  times  our  imple¬ 
mentation  and  TriTri  agreed,  either  by  saying  both  triangles  interseeted  or  both  did 
not,  and  eounted  the  number  of  times  the  2  disagreed. 

Out  of  the  1  million  pairs  of  triangles,  our  GJK  implementation  agreed  with  TriTri 
for  all  but  2  pairs  of  triangles.  In  both  oases,  our  GJK  implementation  reported 
that  the  triangles  interseeted  whereas  TriTri  reported  they  did  not.  Upon  further 
investigating,  we  found  that  in  both  oases  the  triangles  interseeted  at  exaetly  1  point 
along  an  edge.  The  interseotion  point  was  not  any  of  the  vertiees,  nor  were  the 
triangles  ooplanar.  We  believe  that  this  might  aotually  be  a  previously  unknown 
bug  in  TriTri  sinoe  TriTri  reports  that  a  pair  of  triangles  interseot  if  they  share  a 
eommon  vertex.  Thus,  TriTri  olearly  oonsiders  2  triangles  to  be  interseoting  if  they 
share  exaetly  1  point.  Yet,  for  these  2  oases,  TriTri  oould  not  deteot  the  interseotion. 
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This  provides  us  with  confidence  that  our  implementation  reports  correct  results. 
However,  we  cannot  fully  declare  that  our  implementation  to  be  correct.  There 
could  be  an  unknown  bug  lurking  in  either  the  algorithm  or  the  implementation. 
As  mentioned  in  the  introduction,  there  are  numerous  difficulties  in  implementing 
GJK  and  implementations  often  disagree  about  small  details. 

5.  Conclusion 

We  have  presented  a  geometrical  approach  to  the  GJK^  algorithm,  which  may  be 
used  to  detect  if  2  arbitrary  convex  objects  intersect.  The  geometrical  approach  is 
generally  simpler  to  implement  and  may  be  more  intuitive  for  some  readers.  How¬ 
ever,  implementing  the  geometrical  approach  requires  one  to  be  extremely  careful 
to  make  sure  assumptions  are  not  violated.  When  we  first  attempted  to  implement 
GJK,  we  experienced  several  difficulties  because  certain  assumptions  were  never 
explicitly  declared.  Part  of  our  purpose  in  writing  this  report  was  to  lay  out  the 
various  assumptions,  why  they  are  needed,  and  how  to  avoid  common  problems. 

We  compared  our  implementation  of  GJK  against  a  program  that  computes  the 
intersection  of  triangles  and  found  good  agreement,  with  disagreement  in  only  2 
out  of  1  million  cases.  In  both  cases,  our  implementation  was  determined  to  operate 
correctly. 

We  are  currently  using  our  implementation  of  GJK  in  our  system  effectiveness  mod¬ 
els.  Within  these  models,  we  need  to  determine  if  2  objects  collide,  where  the  ob¬ 
jects  could  represent  a  vehicle,  a  threat  against  the  vehicle,  or  a  counter-munition 
designed  to  mitigate  the  threat,  etc.  As  such,  these  objects  have  a  wide  variety  of 
shapes  and  orientations.  We  have  found  GJK  to  be  extremely  fast  and  accurate  in 
determining  collisions  among  these  objects. 
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